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Abstract. New observational constraints on the cosmic matter density Dm and an effectively redshift-independent 
equation of state parameter Wx of the dark energy are obtained while simultaneously testing the strong and null 
energy conditions of general relativity on macroscopic scales. The combination of REFLEX X-ray cluster and type- 
la supernova data shows that for a flat Universe the strong energy condition might presently be violated whereas 
the null energy condition seems to be fulfilled. This provides another observational argument for the present 
accelerated cosmic expansion and the absence of exotic physical phenomena related to a broken null energy 
condition. The marginalization of the likelihood distributions is performed in a manner to include a large fraction 
of the recently discussed possible systematic errors involved in the application of X-ray clusters as cosmological 
probes. This yields for a flat Universe, Dm = 0.291)1 q j 2 and Wx = —0.951q’35 (Ict errors without cosmic variance). 
The scatter in the different analyses indicates a quite robust result around Wx = —1, leaving little room for 
the introduction of new energy components described by quintessence-like models or phantom energy. The most 
natural interpretation of the data is a positive cosmological constant with Wx = — 1 or something like it. 
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1. Introduction 

Measurements of the cosmic microwave background 
(CMB) temperature anisotropies (e.g. Stompor et al. 2001, 
Netterfield et al. 2002, Pryke et al. 2002, Scott 2002, 
Sievers et al. 2002, Spergel et al. 2003), the redshift- 
distance relation of type-la supernovae (Riess et al. 1998, 
Perlmutter et al. 1999), the counts of galaxy clusters 
(e.g. Bahcall & Fan 1998, Borgani et al. 2001, Reiprich 
& Bdhringer 2002), etc., suggest that we live in a dark- 
energy dominated Universe during a phase of accelerated 
cosmic expansion. Perhaps the simplest resolution is to re¬ 
sort to Einstein’s cosmological constant A. Postulating a 
constant leaves many questions unanswered, however, re¬ 
lating to the nature of the particle physics vacuum and the 
(approximate) coincidence in the energy density of dark 
energy and dark matter today. For this finely-tuned con¬ 
stant the answer would seem to lie in the initial conditions. 
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An alternative hypothesis is to consider a time-evolving 
dark energy, while assuming that any cosmological term 
is either zero or negligible. For a time-evolving inhomo¬ 
geneous field (see e.g. Ratra & Peebles 1988, Wetterich 
1988, Caldwell et al. 1998, Caldwell 2002) the aim is to 
understand the coincidence in terms of dynamics. 

A central role in these studies is assumed by the phe¬ 
nomenological ratio w = p/pc^ between the pressure p of 
the unknown energy component and its rest energy den¬ 
sity p. In most investigations the parameter space of w 
is restricted to ru > —1 (exceptions are Caldwell 2002, 
Hannestad & Mortsell 2002, and Melchiorri et al. 2002) by 
assuming that the so-called null energy condition of gen¬ 
eral relativity should be fulfilled on macroscopic scales. 
However, for energy conditions no strict mathematical 
proofs exist hitherto, and their validity is not more than 
a conjecture. Therefore, the present investigation takes 
one step back, using w itself to test the energy condition 
and finding the w value which represents the observational 
data best. 
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Fig. 1. The null energy condition (NEC) and the strong en¬ 
ergy condition (SEC) for a flat FRW spacetime at redshift 
z = 0 with negligible contributions from relativistic particles 
in the parameter space of the normalized cosmic density 
of baryonic and non-baryonic matter and the equation of state 
parameter Wx of a presently unknown energy component. The 
NEC and SEC curves are computed with Eqs. 0 and 0, re¬ 
spectively. The sector between the two curves gives acceler¬ 
ated growing scale factors where NEC is fulhlled but SEC is 
violated. In the sector below the NEC curve all energy con¬ 
ditions of general relativity are violated and the scale factor 
shows a super-accelerated increase. In the sector above SEC 
all energy conditions are fulfilled (especially also SEC) and 
no accelerated cosmic expansion is expected. The vertical line 
at Wx = —1 mark models with Einstein’s cosmological con¬ 
stant (A > 0, especially Cold Dark Matter models ACDM) 
and devides the parameter space into the quintessence-like 
sector (—1 < uix < 0) based on dark energy and the super¬ 
quintessence sector (wx < — 1) based on phantom energy. For 
0 < Wx < 1 ordinary energy might be expected in the form of 
Cold Dark Matter (CDM), Hot Dark Matter (HDM) etc. 

Before we describe the test and its application to as¬ 
tronomical data we briefly review the significance of the 
energy conditions in cosmology relevant for the present 
work and their relation to w (Sect. [21). The test is outlined 
in Sect .(21 In Sect. 0the observational material used for the 
test is described. The results are presented in Sect. Eland 
discussed in Sect. El In the following we define the Hubble 
constant in units of h = i/o/(100kms“^ Mpc“^). 

2. Energy conditions and cosmology 

Assumptions on energy conditions form the basis for 
the well-known singularity theorems (Hawking & Ellis 
1973), censorship theorems (e.g. Friedman et al. 1993) 
and no-hair theorems (e.g. Mayo & Bekenstein 1996). 
Quantized fields violate all local point-wise energy con¬ 
ditions (Epstein et al. 1965). In the present investigation 
we are, however, concerned with observational studies on 
macroscopic scales relevant for cosmology where p and p 


are expected to behave classically. Normal matter in the 
form of baryons and non-baryons, or relativistic particles 
like photons and neutrinos satisfy all standard energy con¬ 
ditions. The two energy conditions discussed below are 
given in a simplified form (for more details see Wald 1984 
and Barcelo & Visser 2000). 

The strong energy condition (SEC): p + Sp/c^ > 0 
and p + p!(? > 0, derived from the more general con¬ 
dition > 0, where is the Ricci tensor for 

the geometry and a timelike vector. The simplified 
condition is valid for diagonalizable energy-momentum 
tensors which describe all observed fields with non-zero 
rest mass and all zero rest mass fields except some spe¬ 
cial cases (see Hawking & Ellis 1973). The SEC ensures 
that gravity is always attractive. Certain singularity the¬ 
orems (e.g.. Hawking & Penrose 1970) relevant for prov¬ 
ing the existence of an initial singularity in the Universe 
need an attracting gravitational force and thus assume 
SEC. Violations of this condition as discussed in Visser 
(1997) allows phenomena like inflationary processes ex¬ 
pected to take place in the very early Universe or a mod¬ 
erate late-time accelerated cosmic expansion as suggested 
by the combination of recent astronomical observations 
(see Sect.0|. Likewise, phenomena related to A > 0 and 
an effective version of A whose energy and spatial distri¬ 
bution evolve with time {quintessence: Ratra & Peebles 
1988, Wetterich 1988, Caldwell et al. 1998 etc.) are al¬ 
lowed consequences of the breaking of SEC - but not a 
prediction. However, a failure of SEC seems to have no 
severe consequences because the theoretical description of 
the relevant physical processes can still be provided in a 
canonical manner. Phenomenologically, violation of SEC 
means w < —1/3 for a single energy component with den¬ 
sity p > 0. For w > —1/3, SEC is not violated and we 
have a decelerated cosmic expansion. 

The null energy condition (NEC): p+p/c^ > 0, de¬ 
rived from the more general condition > 0, 

where Gfju is the geometry-dependent Einstein tensor and 

a null vector (energy-momentum tensors as for SEC). 
Violations of this condition are recently studied theoreti¬ 
cally in the context of macroscopic traversable wormholes 
(see averaged NEC: Flanagan & Wald 1996, Barcelo & 
Visser 2000) and the holographic principle (see covariant 
entropy bound, Mclnnes 2002). The breaking of this crite¬ 
rion in a finite local region would have subtle consequences 
like the possibility for the creation of “time machines” 
(e.g. Morris, Thorne & Yurtsever 1988). Violating the en¬ 
ergy condition in the cosmological case is not as danger¬ 
ous (no threat to causality, no need to involve chronol¬ 
ogy protection, etc.), since one cannot isolate a chunk of 
the energy to power such exotic objects. Nevertheless, vi¬ 
olation of NEC on cosmological scales could excite phe¬ 
nomena like super-acceleration of the cosmic scale factor 
(Caldwell 2002). Theoretically, violation of NEC would 
have profound consequences not only for cosmology be¬ 
cause all point-wise energy conditions would be broken. 
It cannot be achieved with a canonical Lagrangian and 
Einstein gravity. Phenomenologically, violation of NEC 
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Fig. 2. The null energy condition (NEC) and the strong en¬ 
ergy condition (SEC) for a flat ERW spacetime as a function of 
redshift 2 and Wx with negligible contribntions from relativis¬ 
tic particles computed with Eq. dJ. Short dashed curves are 
computed with flm = 0.2, continuous curves with flm = 0.3, 
and dashed-dotted curves with flm = 0.4. Above the curves 
the respective energy conditions are fulfilled. 


means te < — 1 for a single energy component with p > 0. 
The sort of energy related to this state of a Friedmann- 
Robertson-Walker (FRW) spacetime is dubbed phantom 
energy and is described by super-quintessence models 
(Caldwell 2002, see also Chiba et al. Yamaguchi 2000). 
For ic > — 1 NEC is not violated, and the sort of energy 
is termed dark energy and is described by quintessence or 
super-quintessence models. 

3. A test of energy conditions on large scales 

The lack of verification of the energy conditions on macro¬ 
scopic scales suggest that we should first test observation- 
ally the degree to which we can presently trust NEC and 
SEC. It will be seen that the test automatically yields 
those values of the present cosmic matter density and the 
equation-of-state parameter of the dark energy which de¬ 
scribe the observational data best. 

Our starting point is a FRW spacetime filled with a 
positive energy density of unknown nature which can be 
described as a perfect fluid (px ^ 0, px > 0, Wx), with con¬ 
tributions from a pressureless non-baryonic and baryonic 
fluid (pm = 0, Pm > 0, Wm = 0), and from relativistic par¬ 
ticles like photons or neutrinos (pr > 0, pr > 0, Wr = 1/3). 

The different fluids ‘know about each other’ through 
their common gravitational effects and through possible 
explicit couplings of one fluid to the others (a phenomeno¬ 
logical treatment of a local energy transfer between two 
cosmic fluids can be found in Gromov et al. 2002). In 
the late Universe, a coupling between ordinary matter 
and the unknown energy component might remain (e.g. 
Amendola 2000), but the resulting effects are difficult to 
distinguish from predictions of general relativity plus a 


cosmological constant (Torres 2002). We thus regard the 
three fluids mentioned above as effectively independent 
substances over the redshift range covered by the astro¬ 
nomical objects used in our tests. Consequently, for each 
of the cosmic substances a local energy balance holds (e.g. 
Rindler 2001), 


- + 3{1 + w)H = 0, (1) 

P 

derived from the twice-contracted Bianchi identity and 
Einstein’s field equations for a perfect fluid. In iQJ, i? is 
the Hubble parameter and a dot denotes a derivative with 
respect to the cosmic time. 

We are aware that these pre-assumptions are already 
quite specific compared to the usually targeted generality. 
However, our main focus is to learn more about the un¬ 
known energy component ‘x’ and the cosmic phenomena 
related to it. In this sense we proceed further and add all 
the energy and pressure sources to get the net equation 
for the pressumed multi-component effective cosmic fluid. 
In this case, the NEC with its constraint on the passive 
gravitational mass density {p-\-p/c^ > 0) reads 


Px + ^ + Pm + Pr + ^ ^ 0, (2) 

or with the equation of state parameter of the unknown 
fluid defined by Px = and Pr = ^PrC^, 

Px + WxPx + Pm + Pr + gPr > 0 . (3) 


Using the normalized energy density of matter (Um), of 
relativistic particles (Ur) and of the presently unknown 
energy component (Ux), Eq. (0 can be recast into the in¬ 
equality 


Wx > 


Um T -f Ux 


( 4 ) 


NEC for Um, Ur > 0, Ux > 0. Similarily, the SEC with 
its additional constraint on the active gravitational mass 
density (p -I- 3p/c^ > 0) corresponds to 


Wx > - 


Uirr -f 2 Ur -f Ux 

3U^ 


( 5 ) 


SEC for Um, Ur > 0 Ux > 0. Note that the inclusion 
of an unknown energy component in the form of a per¬ 
fect fluid naturally extends earlier definitions of the SEC 
where cases with A-like energies were explicitly excluded. 
However, the present definition still follows the basic idea 
of the SEC because its validity guarantees that the ac¬ 
tive gravitational mass of a multi-component cosmic fluid 
always leads to an attractive gravitational effect. 


d 
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where a = z) is the scale factor of the spacetime, 

2 the cosmological redshift, and where Ac^/3 is replaced 
in the Friedmann-Lemaitre equation by the general term 
87rGpx/3. 
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For the more restrictive case of a spatially flat FRW 
geometry and a negligible contribution from relativistic 
particles (current estimates range from = 0.001 to max¬ 
imal 0.05, see Turner 2002) we have 


Wx > 


1 

1 - ’ 


(7) 


NEC for Om > 0 and 0^ -I- Ox = 1, and 


""" - 3(1 - ’ 


( 8 ) 


SEC for Om > 0 and Om + Ox = 1. Only for the unrealistic 
limit Oni —> 0 the NEC restriction o converges to the 
frequently adopted threshold Wx(min) = — 1. For Ox ^ 0 
and thus 0^ —> 1 one gets r(;x(min) ^ — oo and both 
NEC and SEC are always fulfilled. The link between Wx 
and observable quantities like distances, volumes etc. is 
given by equations of the form 

2 

= Ox /(2:, Wx) + Om (1 + z)^ -|- 

Or (1 + z)^ -|- ( 1 — Ox — Off! — Or ) ( 1 -|- Z )^ , (9) 

which relate H{z) to Wx, 0^, O^ and Ox. The redshift- 
dependency of the latter energy component in m can be 
obtained from the integration of o and is 

f{z,Wx) = expja J [1 -b u;x(z')] t^ln(l-b z')| ■ (10) 

For simplicity we concentrate on a redshift-independent 
constant Wx, so that m leads to f{z,Wx) = (1 -b 
z)3(i+“’*). For Wx = —1 we have / = 1 and Ox corresponds 
to the normalized cosmological constant Oa (Fig.^. 

For constant Wx, the redshift-dependency of 0IH1) can 
be found by replacing the present matter density Om by 
Om(-z) = Oi„/[Oi„ + (1 — Oi„)(l -b leading to the 

condition 

(11) 


1 + z > 


- 1 


7 


1-On 

Om 


H{z) 

Ho 


with 7=1 and 1 for NEC and SEC, respectively, and 0 < 
Om < 1 and Wx < — 7 . Above the redshift limit given in 
Eq. (tnil the respective energy condition is fulfilled. Fig.|^ 
suggests that for z 1 (excluding cosmic epochs with 
e.g. dominating scalar fields like in a possible inflationary 
phase) and for reasonable present values, NEC and 
SEC are always fulfilled - independent of the value of Wx- 
It is also seen that the frequently adopted lower NEC- 
threshold of r(;x(niin) = — 1 is recovered at z <C 0, that 
is, in the distant future. Note, however, that at extreme 
redshifts the constancy of Wx is expected to be a poor 
approximation, suggesting a failure of our considerations 
in these extreme z-ranges. 

This simple discussion shows that for the evaluation of 
NEC and SEC in the given restrictive sense mil) one has 
to measure flm and Wx (most importantly at the present 
epoch). In order to get robust observational constraints we 


utilize the complementarity of two approaches which have 
a different sensitivity on Wx and Recent supernova 
(SN) data constrain Wx and the cosmic matter density flm, 
but the results are highly degenerated (Garnavich et al. 
1998, Perlmutter et al. 1999). The degeneracy can be bro¬ 
ken by using abundance measurements of a large sample 
of nearby X-ray clusters of galaxies which give a precise 
estimate of flm (Schuecker et al. 2003) almost indepen¬ 
dent of Wx (see Sect. 14.111 . The results are compared with 
estimates obtained by a complementary approach which 
combines SN, CMB and other data (Sect.^. 


4. Observational data 

4.1. The X-ray cluster sample 

The ROSAT ESC Flux-Limited X-ray (REFLEX) sam¬ 
ple used for the present investigations consists of the 452 
X-ray brightest southern clusters of galaxies with red- 
shifts mainly below z = 0.3. They are extracted with a 
well-known selection function from the ROSAT All-Sky 
survey (Voges et al. 1999) and confirmed by extensive 
optical follow-up observations within a large ESO Key 
Programme (Bohringer et al. 1998, Guzzo et al. 1999). The 
clusters are located in an area of 4.24 sr in the southern 
hemisphere with Declination < 2.5 deg, excluding galactic 
latitudes |5| < 20 deg and some additional crowded fields 
like the Magellanic Glouds (Bohringer et al. 2001). The 
sample is expected to be at least 90% complete. With this 
sample the cluster X-ray luminosity function (Bohringer 
et al. 2002), the spatial cluster-cluster correlation func¬ 
tion (Gollins et al. 2000), its power spectrum (Schuecker 
et al. 2001), and the cosmic matter density (Schuecker et 
al. 2002, 2003) have been determined with unprecedented 
precision. The 426 REFLEX clusters used for the present 
abundance measurements have at least 10 X-ray source 
counts detected in the ROSAT energy band 0.5-2.0keV, 
X-ray luminosities Lx > 2.5 x 10^^ /i“^ergs“^ and X-ray 
fluxes Sx > 3.0 X lO”^^ ergs”^ cm“^ in the energy band 
0.1-2.4keV. 

The values of the cosmological parameters are esti¬ 
mated by comparing the observed redshift histogram of 
the clusters with model predictions. Note that the vari¬ 
ances of the cosmological parameters shown here are larger 
compared to those obtained with the combined analysis 
of both the redshift histogram and the fluctuation power 
spectrum as used in Schuecker et al. (2003). A brief de¬ 
scription of the model fits of the redshift histograms is 
given below. More details especially the values of impor¬ 
tant model parameters can be found in Schuecker et al. 
(2002, 2003). Deviations from the assumed values intro¬ 
duce possible systematic errors leading to a comparatively 
large scatter of, e.g., the measured as values (see Pierpaoli 
et al. 2002), where as gives the standard deviation of the 
matter density fluctuations in spheres with a comoving 
radius of 8h~^ Mpc. In the present investigation we take 
these systematic errors into account by marginalizing the 
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Fig. 3. Likelihood contours (1-3 (t levels for two degrees of freedom) for various cosmological parameters obtained with the 
abundances of the REFLEX clusters. Note that the lower right panel includes marginalization over as = [0.70,0.95]. The 
parameter priors of each diagram and the marginalization range are given in the main text. 


cluster likelihood distributions over a large as range (see 
below). 

The number of clusters expected under REFLEX con¬ 
ditions at a specihc redshift is given by the integral over 
the mass function where the lower mass limit is a function 
of redshift, flux-limit, cosmology etc. The mass limit is di¬ 
rectly related to an X-ray luminosity via the mass/X-ray 
luminosity relation of galaxy clusters. We use an empir¬ 
ical estimate given in Reiprich & Bohringer (2002). The 
resulting mass limit is transformed into the mass system 
defined in Jenkins et al. (2001) so that the correspond¬ 
ing mass function can be integrated to give the expected 
average number of clusters. 


The computation of the X-ray luminosities takes into 
account the systematic underestimation of the observed 
unabsorbed X-ray fluxes of REFLEX clusters relative to 
the total fluxes (10%, see Bohringer et al. 2002) and the 
cosmic iL-corrections obtained from a rehned Raymond- 
Smith code (Bohringer et al. 2000). For the transforma¬ 
tion of the cluster masses defined at different overdensity 
radii the Navarro et al. (1997) mass density prohle is used 
with a redshift and mass-independent concentration pa¬ 
rameter of c = 5. Deviations from this value in the range 
4 < c < 6 have effects below a few percent and are ne¬ 
glected (Schuecker et al. 2003). We thus assume that the 
REFLEX clusters do not show any signihcant evolution up 
to z = 0.3 as suggested by the redshift-independent distri- 
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bution of the comoving number densities of the REFLEX 
clusters (Schuecker et al. 2001, see also the discussion in 
Rosati et al. 2002). The integration of the mass function 
includes a convolution which takes into account the intrin¬ 
sic scatter of the mass/X-ray luminosity relation, and the 
random flux (luminosity) errors of the REFLEX clusters 
as given in Bohringer et al. (2001). In the present inves¬ 
tigation we use an effective scatter with a formal value 
of 25 % which includes the contributions from flux errors 
(10%) and intrinsic scatter (20%). The intrinsic scatter is 
in reasonable agreement with the Reiprich & Bohringer 
(2002) mass/X-ray luminosity relation if one takes into 
account that realistic mass errors are expected to be a 
factor 1.5 larger than the formal errors given in Reiprich 
& Bohringer (see Schuecker et al. 2003 for a more detailed 
discussion). 

The computation of the theoretical mass function as¬ 
sumes a matter power spectrum and a model for the crit¬ 
ical density contrast which defines the virial cluster mass. 
For the matter power spectrum a Cold Dark Matter trans¬ 
fer function with a given contribution of baryons is as¬ 
sumed (see Eisenstein & Hu 1998). No Wx-dependent cor¬ 
rections of the transfer functions as given in Ma et al. 
(1999) were applied because of the comparatively small 
redshift range (z < 0.3) and scale range (< 1.5 h~^ Gpc) 
covered by the given X-ray cluster sample. For the deter¬ 
mination of the Wx (and Stm) dependent critical density 
contrasts we follow the formalism of Wang & Steinhardt 
(1998) and integrate numerically the relevant ordinary 
differential equations describing the collapse of a spher¬ 
ical over density in an expanding Universe with Wx ^ —1. 
We found that the C function introduced in Wang & 
Steinhardt (their equation All) used to compute the av¬ 
erage critical density contrast between cluster and back¬ 
ground should be replaced by 

c = 5.4[H^(zta)]-°-®®®+0.11[H^(zta)]-'-®''(l-u;x)-'-' , (12) 

giving a good approximation (better than 5%) for the 
range —5 < Wx < 0 at turn-around (^ta)- Eg. (IT^ de¬ 
viates from the corresponding equation obtained for the 
smaller range — 1 < Wx < 0 given in Wang & Steinhardt 
by about 5% and is larger by a factor 1.3 at Wx = —5 (for 
Uni = 0.3). For completeness we also give the equivalent 
linear overdensity of a virialized spherical shell, 

(5c = 1.686^ (^3) 

which is needed when the Press-Schechter and the Sheth- 
Tormen mass functions are used. In the present case 
we use the Jenkins et al. mass function because of its 
higher precision where Jc is not needed. Nevertheless, for 
—5 < ic < — 1 there is almost no change in be except for 
ridiculously small values of Um. 

The comparison of expected and observed cluster 
abundances assumes a Gaussian likelihood distribution. 
For large sizes of the count cells (as used here) this as¬ 
sumption is justified at a high level of statistical signif¬ 
icance by the measurements of Schuecker et al. (2002, 
2003). 


The likelihood contours shown in Fig.O illustrate the 
sensitivity of the REFLEX sample to specific cosmological 
parameters for a flat cosmic geometry. The set of default 
parameter values are: h = 0.70 (Freedman et al. 2001), 
(Tg = 0.711 (normalization of the matter power spectrum, 
see Schuecker et al. 2003, Allen et al. 2002), Um = 0.341 
(see Schuecker et al. 2003), ns = 1-0 (spectral index of 
initial scalar fluctuations, see the recent GMB measure¬ 
ments given in Sect.^l, = 0.022 (baryon density, 

see GMB), Wx = —1. Each panel in Fig.|31shows the l-3a 
likelihood contours for two parameters whereby the re¬ 
maining parameter values are fixed by the default values 
given above. 

Notice that the cosmological parameters erg and espe¬ 
cially Um can be obtained with nearby cluster samples al¬ 
most independent of the value of Wx, as seen in the lower 
panels of Fig.|21 Um measurements based on the abun¬ 
dance of nearby clusters thus appear quite stable against 
the presence of ordinary or more exotic A-like energies. 

The lower right panel of Fig.|21gives the final (wx, Um) 
likelihood values which will be used for the combination 
with the SN data (see Sect. 0. The likelihood distribu¬ 
tion is obtained with h = 0.70, ns = 1, = 0.022 

and after marginalization over the quite large erg range 
of [0.70,0.95]. This interval covers most of the erg values 
obtained recently with different samples, model assump¬ 
tions, and methods (weak leasing, optical clusters. X-ray 
cluster temperature and luminosity functions, power spec¬ 
trum, Sunyaev-Zel’dovich effect power spectrum, galaxy 
clustering etc.) as summarized in, e.g., Pierpaoli et al. 
(2002). Compared to this marginalization range, the rel¬ 
atively small systematics introduced by the uncertainties 
in the values of ns, h, and Ubft-^ are neglected (see Fig.|31 
see also Tablet in Schuecker et al. 2003). 

4.2. The type-la supernova samples 

The two SNela samples described in Riess et al. (1998) 
and Perlmutter et al. (1999) are used for the present in¬ 
vestigations. The distance moduli of the SNe are deter¬ 
mined with spectral and photometric observations. The 
host-galaxy subtracted SN peak magnitudes are corrected 
by the two teams for the cosmic A'-effect, absorption in 
the Galaxy and host galaxy assuming a standard galac¬ 
tic absorption law and no colour evolution (Perlmutter et 
al. did not correct for absorption in the host galaxy but 
made some additional checks and rejected obviously red¬ 
dened SNe from the cosmological tests), time dilation in 
the light curve, and the shape of the light curve. There is a 
distinct difference in the treatment of the latter correction 
in the two samples which complicates a direct comparison 
between the two data sets. The difference is, however, not 
relevant for cosmological applications. 

Gorrections of the change in the peak absolute lumi¬ 
nosity are performed for each SNIa in the Riess et al. 
sample with the Amis method of Phillips (1993), Hamuy 
et al. (1995), and Phillips et al. (1999), and with the Multi- 


Schuecker et al.: Observational Constraints on General Relativistic Energy Conditions 


7 





Wx 


Fig. 4. Likelihood contours (l-Scr levels for two degrees of freedom) obtained with SNela only (left panels) and with SNela plus 
REFLEX X-ray clusters (right panels) marginalized over as, — [0.70, 0.95]. The SN results of the hrst two rows are based on the 
Riess et al. (1998) sample (MCLS corrections upper row, Amis corrections middle row), whereas the data in the lower row are 
based on the Perlmutter et al. (1999) sample. Vertical dashed lines at Wx = —1 represent models with Einstein’s cosmological 
constant. Curved dashed lines devide the parameter spaces into sectors where SEC is valid indicating a decelerated cosmic 
expansion (above the upper dashed curves), sectors where SEC is violated but NEC is fulhlled indicating an accelerated cosmic 
expansion (between the dashed curves), and sectors where NEC (and thus SEC) are violated indicating a super-accelerated 
cosmic expansion (below the lower dashed curves). 


Color Light Curve Shape (MCLS) method of Riess et al. 
(1998). A simple re-normalization of the observed appar¬ 
ent peak magnitude is performed for each SNIa in the 
Perlmutter et al. sample using the stretch factor intro¬ 
duced in Perlmutter et al. (1995, 1997). 

The Riess et al. sample consists of 27 nearby SNela 
{z < 0.2), 10 High -0 SNela (0.30 < z < 0.97) and 6 
High-z Snapshot SNela (0.16 < z < 0.83). The resulting 
catalogue used here gives the redshifts and the distance 


moduli as obtained with the Amis and with the MCLS 
methods. The magnitudes are corrected for the various ef¬ 
fects mentioned above. The errors of the distance moduli 
and cosmological redshifts of the SNe are used for weight¬ 
ing. 

The Perlmutter et al. sample consists of 38 SCP SNe la 
(0.172 < z < 0.830) and 16 Calan/Tololo SNela (0.014 < 
z < 0.101). The sample is the same as used by Perlmutter 
et al. for their Primary Fit (C-fit). Note that this sam- 
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pie excludes 6 SNe from the original sample of 60 SNe 
(2 residual outliers, 2 stretch outliers, 2 likely reddened 
SNe). The resulting catalogue used here gives the SN red- 
shifts, the effective (corrected) peak apparent magnitudes 
and the total uncertainty of the magnitudes. The latter 
quantities already include the uncertainties related to the 
expected errors of the cosmological redshifts of the SNe. 

Note that the 16 low redshift Calan/Tololo SNe la 
and two distant SNe are members of both SN samples. 
Therefore, we cannot regard the Riess et al. and the 
Perlmutter et al. SN sample as statistically independent. 

The present SN likelihood analyses assume Gaussian 
likelihood functions as verified in Riess et al. and 
Perlmutter et al. With the corrected peak magnitudes and 
redshifts of the SNe it is straightforward to compare ob¬ 
served and model magnitudes assuming different values 
of Wx and (including marginalization over method- 
specific quantities as in Riess et al. and Perlmutter et al.). 
The results are shown for the different samples and light 
curve corrections in the left panels of Fig. 01 

5. Combined constraints from X-ray clusters and 
SNe-la 

The combination of the likelihood distributions obtained 
with the X-ray clusters and SNe makes the realistic as¬ 
sumption that both samples are statistically independent 
so that the cluster and SN likelihoods can be point-wise 
multiplicated. The resulting joint likelihood distributions 
of the constraints on Wx and obtained for the three SN 
samples combined with the REFLEX results are shown in 
the right panels of Fig.0 (see also Table 1). As discussed 
in Sect .01 vertical dashed lines at Wx = — 1 represent the 
case of a cosmological constant and devide the parameter 
spaces into the dark energy sectors with — 1 < Wx < 0 and 
the phantom energy sectors with Wx < —1. The curved 
dashed lines are computed with Eqs. 0 and (0 and give 
the deviding lines for NEC (lower dashed curves) and SEC 
(upper dashed curves). 

Figure01 shows that for all three combinations of X- 
ray cluster data obtained with the marginalization inter¬ 
val as = [0.70,0.95] and SN data the centroids of the 
joint likelihood distributions range between —1.05 < Wx < 
—0.75 and 0.28 < < 0.30 (see also Table 1). 

Figure01 illustrates the stability of the results 
by showing the combined likelihood distributions 
for the as marginalization intervals [0.70,0.75] (first 
row), [0.70,0.80] (second row), [0.70,0.85] (third row), 
[0.70,0.90] (fourth row). These computations thus illus¬ 
trate the effects of increasing systematic errors without 
identifying the exact sources of the systematics. In all 
cases the likelihood distributions have maxima at Wx val¬ 
ues between —1.10 and —0.75, and flm values between 0.28 
and 0.32. 

The differences seen in Eigs.0]and01are attributed to 
the different methods used to correct the SN peak magni¬ 
tudes, SN sample-to-sample variations, and different ran¬ 
dom plus systematic errors in the reduction of the cluster 


Tablet. Constraints on flm and Wx, and their la errors ob¬ 
tained with SNe plus X-ray clusters of galaxies after marginal¬ 
ization over as = [0.70, 0.95] assuming a flat geometry. SNe: 
supernova sample, LCC: light curve correction, MCLS: multi¬ 
color light curve shape correction, Amis correction, s: stretch 
factor correction. 


SNe 

LCC 

rim 

Wx 

Riess 

Riess 

Perlmutter 

MCLS 

Amis 

s 

0.281H? 

o.soin® 

o.sotni 

' ^-0.35 

i.UO_Q 35 

i.UO_Q 35 


Average 0.29t[];'(| -0.95l]];^^ 


data. In all cases the centroids clearly fall between the 
NEC and SEC lines. Formal averages over the mean val¬ 
ues and their la errors (without cosmic variance) obtained 
with the largest marginalization range of as = [0.70, 0.95] 
give the final (most conservative) results, 

= 0.291H2. w;x = -0.95t°:i^°. (14) 

Due to the statistical dependencies of the individual SN 
samples, m gives the mean errors obtained with the in¬ 
dividual cluster-SN likelihood combinations, and not the 
errors of the averaged flm and Wx values. 

6. Summary and conclusions 

The null energy condition (NEC) and the strong en¬ 
ergy condition (SEC) of general relativity are tested and 
give observational constraints on cosmic phenomena like 
quintessence, super-quintessence, and Einstein’s cosmo¬ 
logical constant. In order to test NEC and SEC on cos¬ 
mic scales we assume a flat geometry and that baryonic 
matter, non-baryonic matter, relativistic matter and A- 
like matter can be regarded as independent perfect fluids. 
The resulting inequalities (Eqs.[7| and |HJ reproduce the 
frequently used threshold rux = only in the unrealistic 
limit of a zero density of ordinary matter. Note that the 
usually adopted threshold is based on the unknown com¬ 
ponent fix only, which is not enough to cover the large 
diversity of phenomena caused by a mixture of different 
cosmic fluids expected to fill the present Universe (Fig.0. 

In this sense simple cosmological tests of NEC and 
SEC are formulated and applied to the presently largest 
homogeneously selected sample of X-ray cluster of galax¬ 
ies (REFLEX). The most important advantage of X-ray 
clusters as summarized in Borgani & Guzzo (2001) is that 
in contrast to optically selected clusters, their selection 
function (e.g. the sample volume, see Bohringer et al. 
2002) and the relation between cluster X-ray luminosity 
and total cluster mass (see Reiprich & Bohringer 2002) is 
well-known without the need of extensive numerical sim¬ 
ulations. Moreover, the relation between the presence of 
a cluster in X-rays and a peak in the underlying cosmic 
mass distribution (cluster biasing) can be derived from 
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Fig. 5. Likelihood contours (1-3 (t levels for two degrees of freedom) obtained with SNela plus REFLEX X-ray clusters illus¬ 
trating the stability of the SEC and NEC tests by changing the as, ranges for marginalization. First column: Riess et al. (MCLS 
corrections) plus REFLEX. Second column: Riess et al. (Amis corrections) plus REFLEX. Third column: Perlmutter et al. plus 
REFLEX. First row: as = [0.70,0.75]. Second row: as = [0.70,0.80], Third row: as = [0.70,0.85], Fourth row: as = [0.70,0.90]. 


first principles (e.g. Kaiser 1994, Matarrese et al. 1997, 
Moscardini et al. 2000). Therefore, the abundances of the 
nearby clusters are quite sensitive to and almost inde¬ 
pendent of tCx, which is optimal to break the degeneracy 
between flm and shown by the SN data. 

Recently, detailed studies showed that for the applica¬ 
tion of X-ray clusters as cosmological probes several sys¬ 
tematic errors could quite strongly affect the final results 


(H. Bdhringer et al., in preparation, Pierpaoli et al. 2002, 
Reiprich & Bdhringer 2002, Schuecker et al. 2003). This is 
reflected in a large scatter of as values published recently 
by different groups using different samples and methods 
and by comparing X-ray cluster results with results ob¬ 
tained with, e.g., weak lensing, optical clusters. X-ray 
cluster temperature and luminosity functions, power spec¬ 
trum, Sunyaev-Zel’dovich effect power spectrum, galaxy 
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clustering etc. - although some kind of convergence to spe¬ 
cific Oin and as values emerges. The present investigation 
takes this imprecise knowledge into account by a marginal¬ 
ization over a quite large range of erg values [0.70, 0.95] 
which includes about 75% of the values obtained within 
the past two years. In comparison to this range, addi¬ 
tional marginalization over our imprecise knowledge of 
cosmological parameters in the ranges 0.64 < h < 0.80, 
0.018 < < 0.026, and 0.8 < ns < 1-2 was analysed 

and shown to be of secondary importance and is thus ne¬ 
glected in the final results presented here. 

It is quite important to note that the final flni and 
Wx values obtained in the present investigation are al¬ 
most independent on the assumed marginalization ranges. 
Different choices in erg yield changes in the final results 
always smaller than about 7% in both Dm and Wx. An in¬ 
creased marginalization range soley increases the final er¬ 
ror bars and leaves the centroid values almost unchanged. 
We attribute the robustness of the test to the complemen¬ 
tarity of SN-Ia and X-ray cluster data. Future investiga¬ 
tions of Wx will clearly benefit from this complementarity. 

The present analysis neglects the Wx-dependency of 
both the shape, and the amplitude growth of the power 
spectrum of the matter density fluctuations. This is justi¬ 
fied by the limited spatial scale and redshift ranges covered 
by the REFLEX cluster sample. Future cluster samples 
will hardly reach 5 h~^ Gpe scales or so where linear the¬ 
ory Wx-dependent effects on the shape of the power spec¬ 
trum are formally expected. However, Wx can change the 
amplitude of the power spectrum by a factor of about two 
between redshift zero and z = 1 which could be measured 
if the cluster X-ray luminosity/mass conversion and the 
effective biasing parameter of the sample could be com¬ 
puted with high enough accuracy. Deep and wide X-ray 
cluster samples could thus use in addition to the mean 
cluster abundance as the traditional cluster criterion (as 
described in Haiman, Mohr & Holder 2001) another quite 
strong Wx-dependent criterion related to the fluctuations 
of the cluster counts around their mean abundance. 

In the present investigation the joint likelihood distri¬ 
butions of the combination of X-ray cluster and SN data 
are computed as a function of the present-day flm and 
Wx values. NEC and SEC are thus effectively tested at 
redshift zero. 

The combined data fall with about 1.5 — 4cr statisti¬ 
cal significance (depending on the SN sample and light 
curve correction) below the SEC threshold. The SEC can 
thus tendatively be regarded as broken on cosmic scales 
at 2 = 0. A similar breaking of SEC was found by Visser 
(1997) under more general assumptions for the time be¬ 
tween the epoch of galaxy formation and the present where 
he compares the relation between the age of the Universe 
and the age of the oldest observed stars. However, our the¬ 
oretical expectations suggest that the broken SEC state 
should not hold for redshifts z > (0.28-1.20) lEa.ITTIand 
the Icr error corridor given in Eq. m), in contrast to the 
larger z range implied by the analysis of Visser. 


The data fall above the NEC threshold with about 
1.5 — 3(7 statistical confidence, again depending on the 
SN sample and light-curve correction used. NEC can thus 
tendatively be regarded as fulfilled. 

In every case tested sofar, the combined X-ray cluster 
and SN data obviously populate the sector between NEC 
validation and SEC violation and thus provide further ob¬ 
servational evidence for an accelerated cosmic expansion 
at z = 0. 

The observational constraints obtained from the com¬ 
bination of X-ray cluster and SN data on Dm are in 
good agreement with recent cluster data (e.g. Borgani et 
al. 2001, Allen et al. 2002, Schuecker et al. 2002, 2003, 
Pierpaoli et al. 2002), constraints from CMB data (see 
references given in Sect.^ especially the WMAP result 
Dm = 0.29 ± 0.07, 68% confidence, WMAP data only, 
Spergel et al. 2003), and galaxy data (Szalay et al. 2001, 
Lahav et al. 2002). 

Our results on Wx are consistent with the constraints 
obtained from type-la supernovae (Carnavich et al. 1998, 
Perlmutter et al. 1999), with recent CMB data (e.g., 
Baccigalupi et al. 2002, Bond et al. 2002) and with 
the baryonic fraction in galaxy clusters (Ettori et al. 
2003). The results obtained by the combined analysis of 
CMB and SN data of Hannestad & Mortsell (2002) and 
Melchiorri et al. (2002) yield the respective 95% confidence 
constraints —2.68 < Wx < —0.78 and —1.62 < Wx < —0.74 
(see also Caldwell 2002), quite consistent but slightly 
larger than the results obtained with the SN and X-ray 
cluster data given here. Finally, the constraint Wx < —0.78 
is obtained with 95% confidence from the combination of 
WMAP, SN, 2dFCRS and Lya data (Spergel et al. 2003). 

Type-la SN and X-ray cluster data thus support a pic¬ 
ture of a universe which is presently in a state of acceler¬ 
ated expansion, where NEC is most probably not violated 
(no super-acceleration with a possible catastrophic end¬ 
ing, MeInnes 2002), and in which a cosmological constant 
or something like it provides the dark energy. 

However, one still has to be cautious with conclusions 
about NEC and SEC because they are only tested un¬ 
der restricted conditions. Furthermore, we are aware of 
the necessity to study in much more detail our assump¬ 
tions that both galaxy clusters and SNe la do not evolve 
over the respective redshift ranges covered by the given 
observations. Moreover, important relations like the total 
cluster mass/X-ray luminosity relation and the extinction 
law relevant for nearby and distant SNe have to be known 
with much higher precision because the deviations from 
Wx = ~1 at z = 0 for interesting dark energy scenarios 
might be smaller than the error bars of the present results. 
The present investigation tries to take into account possi¬ 
ble systematic errors of the treatment of the X-ray clusters 
by using a quite large as marginalization interval. Future 
measurements based on improved relations can work with 
smaller intervals expected to provide quite precise cosmo¬ 
logical constraints on both D^ and Wx- 
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